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<^ ■ We evaluate the local contribution g^L of coherent matter with lagrangian den- 

sity L to the vacuum energy density. Focusing on the case of superconductors obey- 
ing the Ginzburg-Landau equation, we express the relativistic invariant density L in 
terms of low-energy quantities containing the pairs density. We discuss under which 
physical conditions the sign of the local contribution of the collective wave function 
to the vacuum energy density is positive or negative. Effects of this kind can play 
an important role in bringing about local changes in the amplitude of gravitational 
vacuum fluctuations - a phenomenon reminiscent of the Casimir effect in QED. 



X 

j^ ■ 1. Global cosmological term 

The vacuum energy density term, or "cosmological term" in the Einstein equations 
has become popular again in the last years, after a series of observations which indicate 
A ~ 10" 50 cm~ 2 as most probable value. Earlier observations set an upper limit on A of 
the order of 10~ 54 cm~ 2 , so it was thought to vanish exactly for symmetry reasons. 

As remarked by several authors, the observed non-zero value of A creates an arduous 
fine-tuning problem pQ . This value should in fact be regarded as the residual of a complex 
interplay, in poorly known sectors of particle physics, between positive and negative vacuum 
energy densities. According to [2j, the observed residual should also be scale-dependent 
and this dependence could appear most clearly at length scales corresponding to the mass 
of the lightest particles, like neutrinos or unidentified scalars. 

The above mentioned value of A is global, in the sense that it is supposed to define 
a uniform background present in the whole universe. In order to elucidate the notion of 
local contributions to the vacuum energy density, let us first recall here some known facts. 
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The cosmological term is an addition to the l.h.s. of the Einstein equations, of the form 
Ag^x), with A constant: 

R^( x ) - 7,9^( X ) R ( X ) + A 9^( x ) = SnGT^x). (1) 

A quantum field $ contributes to the cosmological constant in two ways: 

(i) Its zero-point oscillations are associated with a huge energy density, whose value 
depends on the frequency cut-off. 

(ii) The field can have a non- vanishing vacuum expectation value (VEV) $o — (0|^|0)- 
This contribution is essentially classical. If L is the field lagrangian, the corresponding cos- 
mological term is — 87iGL(&o)g flu (x), because the energy-momentum tensor has in general 
the form 

T^ = n M d„$ - g^L, (2) 

where IT^ (conjugate momentum of $) and d v Q vanish in the case of a constant field. One 
of the roots of the cosmological constant paradox is clearly seen here: any constant term 
in the lagrangian of a system, though irrelevant for the system's dynamics, has gravita- 
tional consequences. In elementary particle physics a non-vanishing VEV is usually the 
consequence of a spontaneous symmetry breaking process. 

2. Local cosmological term 

A local contribution to the vacuum energy density can arise when the state of a 
localized physical system is described by a classical field comparable with the VEV of a 
quantum field. We are interested into cases of this kind occurring in condensed matter 
physics. In this context, the physical systems properly described by continuous classical- 
like fields (also at microscopic level, not just in a macroscopic-average sense as for fluids) 
are basically 

(1) the electromagnetic field in the low-frequency limit, in states where the photons 
number uncertainty is much larger than the phase uncertainty; 

(2) systems with macroscopic quantum coherence, described by "order parameters", 
like superfluids, superconductors and spin systems. 

Suppose that, in one of these systems, the field has a constant value $o i n a bounded 
region and is zero outside. (Consider for instance a container with superfluid helium of 
constant density.) We can speak of a contribution of the field to the cosmological constant 
in this region, equal to — 8ttGL(Qq). A magnitude order estimate shows that such a 
contribution can be of the same order of the observed A or even larger, although it is 
clear that this energy, being present only in a small region of space, does not have any 
cosmological relevance. 

From the classical point of view, one can correctly object that the description of this 
situation in terms of a local cosmological constant is purely formal, because the gravita- 
tional field present is just that due to the superfluid regarded as an energy- momentum 
source. Moreover, there is no distinction, still at the classical level, between a truly con- 
tinuous source, like the superfluid wave function, and an incoherent fluid. 

3. The cosmological term vs. dipolar fluctuations 

The perspective changes if one takes into account short-scale gravitational quantum 
fluctuations. Suppose to describe gravity with the covariant perturbation theory in the 



weak-field approximation. The action contains some parameters, and one of these is the 
effective A in the considered region. The value of A acquires a special significance in relation 
to the so-called "dipolar vacuum fluctuations" , which constitute a class of zero modes of 
the pure Einstein action [3|. Let us recall their features in short. 

The functional integral of pure Einstein quantum gravity can be written as z = 

f d[g !_,,„] exp(iS/h), with S = J d 4 xJg(x)R(x)/87rG. The "spacetime foam" [I] consists 

of fluctuations whose action does not exceed a quantity of order ft. 2 This implies, for 
curvature fluctuations on a scale r, that \R\ < G/r A (according to naive power counting) or 
\R\ < l/(Lpr) (according to numerical lattice estimates jS]). Therefore, large fluctuations 
are expected to take place only at very small distances. 

Since, however, the Einstein action is not positive definite, one can also expect some 
fluctuations due to peculiar cancellations of distinct contributions to the action, which are 
by themselves larger than h. In order to show this explicitely, let us consider the Einstein 
equations with an auxiliary source T^: 

R^u{x) - -g lu/ (x)R(x) = -8nGT^(x) (3) 

and their trace 

R(x) = 87iGg^{x)T^{x) = 8nGTiT{x) (4) 

Then consider a solution g^ v (x) of (|HJ) with any source satisfying the condition 

/ d i x^g{x)i:iT(x) = (5) 

In view of (jlj), this metric has zero action. We have constructed in this way a zero mode 
of the pure Einstein action. The source is unphysical, but it is "forgotten" after obtaining 
the metric. Condition (J3J) means in fact that it is a "dipolar" source, with a compensation 
between regions having positive and negative mass-energy density. Since this auxiliary 
source is used to construct a virtual field configuration, we sometimes call it a "virtual 
source" . 

In the presence of a cosmological term in the action, the dipolar fluctuations acquire 
a small monopolar component. In fact, the contribution of a dipolar fluctuation to the 
cosmological term is of the form [3 

AS A = ArMr 2 Q, (6) 

where r is the duration of the fluctuation, M is the order of magnitude of the virtual +/- 
masses, r their distance and Q is an adimensional function which can be positive or negative 
and depends on the detailed form of the dipole. This means that it is still possible to obtain 
a zero mode, provided the contribution to the Einstein part of the action is the opposite of 
AS\. This implies in turn that the integral of the covariant trace of the energy-momentum 
tensor needs not be zero, i.e., there is a monopolar mass component in the virtual source. 

The magnitude order of the monopolar component is readily found to be 

|m iot |~AMr 2 (7) 

Here the amplitude M of the dipolar component is limited only by the i? 2 -like quantum 
corrections to the Einstein action (Ref. [3], Section 2.4). Choosing length and time scales 

2 Note, however, that we use units such that h = c = 1 throughout this work. 



as r ~ 1 cm, r ~ 1 s (~ 10 10 cm in natural units), one finds an upper limit for M of ~ 
10 56 cm -1 ; with the observed value of A, the monopolar residual is then \mtot\ ~ 10 6 cm" 1 . 

So we see that the effective local value of the cosmological constant sets the am- 
plitude of the monopolar component of a class of strong quantum gravity fluctuations, 
corresponding to remarkable positive and negative virtual mass densities, and present at 
distance scales much larger than Planck scale. 3 Consider now the local contribution of 
coherent matter to the cosmological term again . If this causes some local inhomogeneity in 
A, an inhomogeneity in the vacuum fluctuations will follow, and it is well known from the 
Casimir effect that this can give rise to observable effects. In an analogy with the Casimir 
effect, the role of the metallic conductors which impose a cut-off on the electromagnetic 
vacuum fluctuations is played here by the coherent condensed matter. 

4. General form of L for a superconductor obeying the Ginzburg-Landau equa- 
tion, in the non-relativistic limit 

At this point it is important to evaluate the lagrangian density L in some concrete 
situation and compare its value, in sign and magnitude, with the background A/8ttG. For 
the electromagnetic field this is straightforward. Superfluids, on the other hand, cannot be 
adequately described by an effective lagrangian theory 0. So we shall focus our attention 
on the case of superconductors, for which a widely used lagrangian phenomenological theory 
is available - the Ginzburg-Landau theory - based on the macroscopic wave function ipGL- 
In the following we shall compute the relativistic lagrangian of a superconductor in terms 
of the non relativistic wave function ipci, which is well known in several physical situations. 

We are therefore seeking the low energy limit of a scalar field theory appropriate for 
superconductors. This is not trivial. Some proposals for the inverse procedure, namely a 
relativistic generalization of Ginzburg-Landau theory, have been previously described in 
the literature [SJ E] . However, in only variations in the phase of ipci are considered, 
and in [7j the tetrad formalism is employed, with a particular gauge fixing. We shall 
follow a somewhat more conservative method. Let us consider the usual time-independent 
non-relativistic limit of the Klein-Gordon wave function, namely 



fe(x) = e mt <\>{x), (8) 

where m is the Cooper pairs mass. By introducing this into the Klein-Gordon lagrangian 
one obtains 

L = -^V^ G (x)Vfe(x). (9) 

This coincides with the part of the Ginzburg-Landau lagrangian containing the partial 
derivatives [8J, provided the wave function is suitably normalized. We set 



*Pgl(*) = Vmil) KG (x) (10) 

and so obtain 

i=-^ V ^WVfe(x). (11) 

2m 

In the following we shall just write ip instead of ipGL and we shall omit the space dependence. 
The normalization above corresponds to |-0(x:) | = p(x), where p 2 is the density of Cooper 
pairs. 



3 In a complete, non-perturbative theory of quantum gravity these fluctuations should of course take part 
in defining the renormalized value of the coupling constants and of quantities like the vacuum polarizability. 



Next we generalize our initial lagrangian, without spoiling its relativistic invariance, 
in such a way to make it coincide with the full Ginzburg-Landau lagrangian in the low 
energy limit. To this end it suffices to add a quadratic and a quartic term. We also 
introduce the minimal coupling with the four-potential A(x), finally obtaining 



L = — —| - iVi) + 2eA^| 2 - aip*ip - -p^ip) 2 



(12) 



2m 



(V^*)(V^)-2iAip*V^ + 2iAipVi/j* + Ae 2 A^*^ - atp*ip - -P(ip*i^) % (13) 



where a and (3 are two arbitrary coefficients, which in the end are identified with the 
Ginzburg-Landau coefficients. This lagrangian density, multiplied by — 8nG, gives the 
local contribution to the vacuum energy density. Now we evaluate it for wave functions 
which satisfy the Ginzburg-Landau equation. 

A crucial point should be stressed here. Remember that this equation is obtained by 
minimizing the spatial integral of L; while doing this, one transforms a term (Vil>*)(Vip) 
into a term -0*V 2, ?/>, integrating by parts and supposing that ip = at the boundary. Here, 
however, we are not interested into the integral of L, but into its local value. Therefore, 
we take the Ginzburg-Landau equation 



1 
2m 



-iVip + 2eA) 2 i) + (a + (3ip*ip) ip = 



(14) 



and multiply it by ip* on the left. In the gauge VA = 0, V commutes with A; after 
isolating the term 4e 2 A 2 ifj*tfj and replacing it into (fT3*|) . we find 



L = 

2m 



(W)(V^) - 2ieA(^*Vij + ipVij*) + ^*V 2 ^ - m(3(ip*i)) 



(15) 



We observe that L is, by construction, real. Therefore the imaginary part of the 
expression above must vanish and we have 



2ieA(ip*V^ + VW>*) = iIm(tfj*V 2 tp). 
In conclusion equation (|T5|) becomes 

1 



L = — 



2m 



*„l.\2 



(V</>*)(W0 + i?e(?/>*V» - m(3(ip*ip) 



(16) 



(17) 



5. Numerical evaluation of L 

Introducing variables p(x) and #(x), such that ^(x) = p(x)e l6 '^ x - ) , we finally obtain 
for the lagrangian density a simple expression, in which the magnetic field does not appear 
explicitly 

1 



2m 



(Vp) 2 + P V 2 p - mf3p 4 



(18) 



In this expression we can restore h, obtaining 



L = 

2m 



h 2 {Vpf + h 2 pV 2 p-mPp 4 



(19) 



Table 1: Magnitude orders of the lagrangian densities L\ and L 2 for a type I superconductor 
(Pb) and for a type II superconductor (YBCO), computed according to eq. (fT9"j) . Also 
listed are the values of the London length A and the coherence length £ at T = jH]and 
the average pairs density n p , computed from A by the relation n p = m e /(2/i e 2 A 2 ). For 
YBCO, the values of £ along the a-b direction and the c direction are given separately, and 
so the corresponding values of L 2 ; n p is the density in the a-b planes, computed with the 
effective mass m* ~ 4.5m e . 





Pb 


YBCO 


A (m) 


3.9- lO -8 


1.4 ■ 10~ 7 (a-b direction) 


f O) 


8.2- 10 -8 


1.6 • 10~ 9 (a-b direction) 
2.4 ■ 10~ 10 (c direction) 


n p (m~ 3 ) 


9.3 ■ 10 27 


3.5 ■ 10 27 


In (Jm- 3 ) 


10 4 


10 6 


L 2 (Jm- 3 ) 


10 4 


10 6 (a-b direction) 
10 8 (c direction) 



We recall [S] that the a and (3 coefficients depend on the absolute temperature T. 
The coefficient (3 is always positive and approximately constant near T c ; a is negative 
for T < T c and behaves like ~ const. (T — T c ) near T c . The ratio between a and /3 is 
given by the relation n p = —a//3, where n p is the average density of pairs in the material. 
Finally, (3 is linked to the value of the Ginzburg-Landau parameter k = A/£ by the relation 
k 2 = m 2 (3 1 '(2/jLoh e 2 ). We also recall that the wave function must satisfy suitable boundary 
conditions at interfaces with vacuum and at junctions with normal conductors. 

It is straightforward to check that the sign of L is positive for two types of configu- 
rations: 

(1) For the constant solutions of eq. (|14J1 in the absence of external field, which implies 
p 2 (x) = n p . The corresponding constant lagrangian density is L\ = \(3n 2 p . 

(2) For regions of the superconductor where pV 2 p is negative and greater, in absolute 
value, than (Vp) 2 . It is straightforward to check that these are regions located around local 
density maximums, or more generally lines and surfaces where the first partial derivatives 
of p are zero and the second derivatives are negative or null. The lagrangian density at 
a maximum is L 2 ~ f^p|p"|- If the maximum is sharp, L 2 can be much larger than L\. 
Configurations of this kind are characteristic of solutions of the Ginzburg-Landau equation 
with strong magnetic flux penetration j^J. 

In all other configurations, L is negative. Some numerical estimates are given in 
Table 1, where the gradients are taken to be of the order of p£ -1 . At local minima or in 
regions with strong gradients, we can suppose that \L\ is of the same magnitude order as 
L 2 . 



6. Conclusions 



We are interested into the value of the lagrangian density L for condensed matter 
systems described by an order parameter, because the value of L - by interfering with 
the vacuum energy term in the Einstein action - affects locally the amplitude of dipolar 
vacuum fluctuations. This amplitude is proportional to \L — A/8ttG\, where A is the 
"natural" background vacuum energy density. At cosmological scale A is known to be 



negative (with our conventions on the metric) and of the order of 10~ 50 cm~ 2 , so that 
|A/87rG| ~ 10" 1 J/m 3 in SI units. However, A is probably scale-dependent, and increases 
at small distances |2]. 

For electric and magnetic fields in the low frequency limit, when the photons number 
uncertainty is much larger than the phase uncertainty, the lagrangian density is simply 
proportional to the square of the field strength and has positive and negative sign, respec- 
tively. The magnitude order of L is, for instance, Le ~ 10 2 J/m 3 for E = 10 6 V/m and 
L B 10 5 J/m 3 for B = 1 T. 

For a superconductor described by a collective wave function obeying the Ginzburg- 
Landau equation, we have seen that in the non-relativistic limit it is possible to express 
L as a function of I^glI- Relevant values of L are those for the case of constant pairs 
density (Li) and those at local density maximums (L2). As shown in Table 1, this kind 
of local contributions to the vacuum energy density can be much larger than the pure 
electromagnetic contributions. 

In the absence, however, of any reliable experimental or theoretical estimates for the 
value of A at short scales, it is not possible to compare A/8nG with the above values of L 
and guess whether the lagrangian density of coherent matter can affect the amplitude of 
the monopolar component of dipolar vacuum fluctuations. Conversely, the experimental 
search for anomalous effects of the gravitational fluctuations in the presence of coherent 
matter can give useful information on the short-scale value of A (see ^1] and ref.s). 

The hypothesized phenomenon could be regarded as a sort of gravitational analogue 
of the Casimir effect. In spite of the relative weakness of gravitational forces at atomic scale, 
such vacuum effects can be relevant because the dipolar fluctuations in Einstein gravitation 
are zero-modes of the action and are therefore much stronger than electromagnetic vacuum 
fluctuations. 
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